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Abstract 



We introduce a time-dependent perturbation method to calculate the number 
of created particles in a ID cavity with an oscillating wall of the frequency Q. 
This method makes it easy to find the dominant part of the solution which 
results from the parametric resonance. The maximal number of particles are 
created at the mode frequency f2/2. Using the Floquet theory, we discuss the 
long-time behavior of the particle creation. 
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I. INTRODUCTION 



The particle creation by the parametric resonance is a very important phenomenon to 
understand the mechanism of reheating after inflation in the early universe [0. The pho- 
ton production by the parametric resonance in a cavity with an oscillating wall is another 
important phenomenon to be observed in the experimental situation. Recently, the photon 
creation in an empty cavity with oscillating boundaries has attracted much attention |fj || . 
It was also proposed that the high-Q electromagnetic cavities may provide a possibility to 
detect the photons produced in the nonstationary Casimir effect ||[7|]. Therein, they con- 
sidered the vibrating wall with the frequency Q = 2oj\ and found the resonance excitation 
of the electromagnetic modes. 

In this paper we consider the quantum (electromagnetic) field in a cavity with an oscillat- 
ing boundary of the frequency Q. We calculate the number of particles (photons) produced 
by the parametric resonance. For the calculation of time-evolution of quantum fields, we 
introduce a time-dependent perturbation method that makes it possible to calculate the 
photon number for any Q and to observe clearly the effect of the parametric resonance. For 
a long-time behavior of the quantum field we use the Floquet theory. 

The organization of this paper is as follows. In Sec. [II] we review the scheme of the field 
quantization in the case of moving boundaries. In Sec. [ITT] we introduce a new perturbation 
method to find the time evolution of the quantum field. Here we write the dominant part 
of the solution of wave equation which results from the parametric resonance. We calculate 
the number of particles created by the vibration of the boundary ||. In Sec. |V| we use 
the Floquet theory in the perturbation expansion to examine the long-time behavior of the 
solutions. We develope the method to find the characteristic exponent of the solution and 
the periodic part of the corresponding solution that is linear combinations of the mode 
functions. We get the three term recurrence relations between the coefficients of the mode 
functions. The last section is devoted to the summary and discussion. 

II. QUANTUM FIELDS IN A ID CAVITY WITH A MOVING BOUNDARY 

Let us consider a quantum field obeying the wave equation (c = 1) 

d 2 A _d 2 A _ 

with time-dependent boundary conditions: 

A(0,t) = = A(L(t),t). (2.2) 
The field operator in the Heisenberg representation A(x, t) can be expanded as 

A(x,t) =£[Mn(M)+teOM)], (2.3) 

n 

where and b n are the creation and the annihilation operators and tp n (x,t) is the corre- 
sponding mode function which satisfies the boundary condition ip n (0,t) = = ip n (L(t),t). 
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For an arbitrary moment of time, following the approach of Refs. PHTTI, we expand the 
mode function as 

lp n (x, t)=J2 Qnk{t)M X > ( 2 - 4 ) 
k 

with the instantaneous basis 



Here Q n k(t) obeys an infinite set of coupled differential equations 

Qnk + ^k 2 {t)Qnk = 2A gkjQnj + ^ X! 9kjQ 



n] 
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where A = L/L and 



+^ 2 J29jk9jiQni (2.6) 

3,1 



9k >-{ o ( 3 = k)- (2 ' 7) 



and the time- dependent mode frequency is 



Mt) = m (2-8) 



For L(t < 0) = Lo, the right hand side of Eq. Q2.6Q vanishes and the solution in this 
region is chosen to be 



e ~iuj k t 



Qnk(t) = ^— 5 nk . ( 2 - 9 ) 

so that the field ( |2.3|) can be written as 

g— iuj n t 

A(x,t < 0) = ]>> n ^=^ n (x,L ) + H.c] (2.10) 

n V ZuJ n 

where u n = j^. For the scalar field or the electromagnetic field, the Hamiltonian has the 
form H = Y^n^niblfin + 5) j an d we can interpret as the number operator associated 
with the particle with the frequency u n . 

After the change of the boundary, we assume L(t > T) = L , then the solution of 
Eq. ( |2.6|) with the initial condition ( |2.9|) , can be written as 

Qnk(t >T) = a nk — + (3 nk — (2.11) 

V V 

From (|2.3|) and (|2.4|), we have 
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e -iu) n t 

A{x, t>T)= ^[a n ^=^„(x, Lo) + H.c], (2.12) 

n V ZuJ n 

where 

a fc = 51[&„a nfc + bi(3* nk ] 

n 

al = YP n al k + b n (3 nk \. (2.13) 

n 

Further, it follows from H = J2n u n{ a n a n + \) that a n a n is the new number operator at 
t>T. 

If we start with a vacuum state |0&) such that b |0&) = 0, the expectation value of the 
new number operator is 

oo 

N k = (0 6 |4a fc |0 6 > = £I/U| 2 , (2.14) 

n=l 

which is the number of created particles. (Note that the quantum state does not evolve in 
time in the Heisenberg picture.) 



III. TIME EVOLUTION OF THE QUANTUM FIELD IN A CAVITY WITH AN 

OSCILLATING BOUNDARY 



In this section we find the time evolution of quantum field operator (|2.3| ) by solving 
Eq. ( |2.6|) with the motion of the wall given by 

L(t) = L [l + esm(Qt)). (3.1) 

Here Q = = '-/tt/Lq and e is a small parameter characterized by the displacement of the 
wall. For e C 1, having in mind that A(i) ~ e and taking the first order of e in the mode 
frequency ( |2.8| ) 

A/7T 

Mt) = -ri 1 + esm(nt)}-\ (3.2) 
we can replace Eq. ( |2.6| ) by a pair of coupled first-order differential equations 

Qnk Pnk 

P nk = -uJ k 2 (l ~ 2e sin Vtt)Q nk + 2y ^ g kj P n j 

L 3 

+TJ29kjQnj + 0(e 2 ). (3.3) 
L 3 

Introducing the new dynamical variables 

X nM = ^(Qnk±i^) (3.4) 
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and the vector notation 



X n (t) 



the above equation (Of) can be written as 



d 



X n ,2- 



V 



(3.5) 



-X n (t) = V^X n (t) + eV {1) X n (t) 



where and are matrices given by 



and 



V ka,jcr> — U ^ V ku,ja' e 



(3.6) 



(3.7) 



(3.8) 



where 



Li — 1. 

k [ 2 +S 4j 



S(T-d kj 



(3.9) 



with s, a, a' = +, — . Here we used Q = ju-i and Uk = kuj\. 

To find the solution of Eq. (|3.6|), we introduce a perturbation expansion: 



X n = Xf + eX« + 6 2 X( 2 ) + 



By inserting ( |3.10| ) into Eq. ( |3 . 6| ) , identifying powers of e yields a series of equations: 

dt n n 

-jfW = v^x^ + v<®x® 

^£ n n 1 r n 



(3.10) 

(3.11) 
(3.12) 



From the initial condition fl2.9|) , we have the solution to zeroth order equation (|3.11|) 

4°L = SrtSo-e***, (3.13) 
and to the first order equation (|3.12|) 



-i(cFk—s^+n)u)\t' 



(3.14) 



When the exponent of exponential function in the integrand of (|3.14j) vanishes, we have terms 
proportional to Uit which are the effects of parametric resonance. In the usual situation, 
since U\t ^> 1, only the resonance terms are dominant and the solution becomes by retaining 
only them: 
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Qnk(t) ~ /pi ^ * $nk 



eujit 

y/2uJk~ l 

k-,n- e lUJkt $k,n+-y + Vk-,n- e ^fc,n--y] ■ 

(3.15) 



After some time interval T the wall stops at x = Lo, then the solution is described by 
( P- By comparing (|2.11 ) with ( |3.15 ), the Bogoliubov coefficient (3 n k can be read from 
the solution Q n k to the leading order in e 

(3 nk = euiTv£ +tn _8k,i-m (3.16) 

which is the coefficient of negative frequency mode function in ( |3.15| ). 
Using (|2.7|) and Q3.9|) , finally we have 

Wnkf = ~nk(eu 1 T) 2 6k,v-n. (3.17) 
Therefore the total number of particles created in the k th mode from the empty cavity is 

n=l ^ 

This result is a generalization of Ref. |7j in the short time limit (ec^T <C 1) and it agrees 
with that result for 7 = 2 and k = 1. It should also be noted that the maximal number of 
photons are created at the mode frequency 

fc=~or ^fc = 77- ( 3 - 19 ) 
for 7 = even and at its nearest neighbor frequencies k — (7 ± l)/2 for 7 = odd. 



IV. LONG-TIME BEHAVIOR OF THE SOLUTION: PERTURBATION 
APPROACH USING FLOQUET THEORY 

In this section we discuss the long-time behavior of the solution to the differential equa- 
tion ( |3.6|) . Although the time-dependent perturbation method developed in the previous 
section gives the method to calculate the higher order solution, it does not provide the con- 
vergency of the solution. So it is difficult to examine the long-time behavior of the solution. 
Here we develop another perturbation method using the Floquet theory. 

Consider a r-periodic system of differential equations such as ( |3.6|) : 

^X(t) = [V^ + eV w }X(t) (4.1) 
etc 

with the periodic condition 

[y(0) + e y(i)](j + r ) = [y(0) + e y«](£), (4.2) 
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where r = 2tt/Q. The Floquet theory states that the solution of (|4.1| ) should be of the 



form [13 1 



X(t) = e^Z(t) 



(4.3) 



where Z(t) is r-periodic. For simplicity of the indices of the matrix, we will slightly change 
the notation: 

( : \ 

X 2 
\ ) 



X(t) 



(4.4) 



V^f = ikui8 k j, 



(4.5) 



and 



s=± 



(4.6) 



where 



v k,j = 19kj' 
and k and j are nonzero integers. 



To find the solution of Eq. (4~1), we introduce a perturbation expansion: 



(4.7) 



x[Z (0) (t) + eZ {l \t) + e 2 Z (2) (t) + ...]. 



(4.8) 



The zeroth order equation 



tit 



and the first order equation 



d Z(i 



dt 



can be easily solved as: 



Zf = C k e iuJkt 



(4.9) 



(4.10) 



(4.11) 



and 
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z (l) = e ik^t J 1 dt ' e -ik^t' 



X 



EE^^-M^e^. (4.12) 



When the exponent of the exponential function in the integrand of (|4.12|) vanishes, the 
integration gives the term proportional to t. This contradicts the periodicity condition of 
Z(t), therefore the coefficient of such term should vanish. Thus, we have the following 
recurrence relation 

^fc+ 7 Cfc+7 - ^ C k + ^fc- 7 Cfc-7 = 0. (4.13) 

Note that the coefficients are coupled to the 7th neighbor modes. This three term recurrence 
relation can be written as the following linear equation 

M(/x 1 )C ? = (4.14) 

From the condition for the existence of the nontrivial solution, we can find the characteristic 
exponents by solving 

detM(/ix) = 0. (4.15) 

Then we have the eigenvalues /i^ and the corresponding eigenvectors C A , where we intro- 
duced the superscript A to distinguish the eigenvectors. Then the characteristic solution 
is 

X£{t > 0) = e ^?^c£e lkwit . (4.16) 

Note that the amplitude (the coefficient of the harmonic function) is exponentially increasing 
when the real part of the characteristic exponent fi\ is positive. This is the effect of the 
parametric resonance. 

By linear combinations of the characteristic solutions, we find the solutions 

X n , k (t) = J2dnAX£(t) (4.17) 

A 

that satisfy the initial conditions 

X n , k {t<0) = e ikuit 6 nk e(-n), (4.18) 

where 9 is the Heaviside unit step function. From the intial conditions, the coefficients d n A 
are obtained by solving the linear equations 

Y< d nAC£ = 5 nk 6(-n). (4.19) 

A 

Since the matrix equation is infinite dimensional, it is difficult to find the full solution in 
a closed form. However, it is helpful to consider the truncated matrix in order to understand 
the long-time behavior of the solution. As the simplest model (7 = 2), we consider only the 
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mode frequency to the next neighboring frequency, e. g. to the for u\ mode, we have the 
following eigenvalue equation: 



( -Hi -y/3/2 
y/3/2 -n± -1 

-1 -/i! y/3/2 

{ -V3/2 -fn t 

then this equation has the following eigenvalues 

H\ = ±(1 ±iV2)/2 
and the corresponding eigenvectors 



\ 












/ 


\ c 3 J 



/ l 

-2/ii/V3 

y/3/2 + 2/4/y/3 

\^ l + A^l-l)/3) ) 



(4.20) 



(4.21) 



(4.22) 



In the long-time behavior, the solutions of which characteristic exponent is positive will be 
dominant. 



V. DISCUSSION 

We developed a perturbation method to find the time-evolution of the field in a cavity 
with an oscillating boundary. This method makes it possible to calculate the particle num- 
ber for any oscillation frequency of the boundary and to observe clearly the effect of the 
parametric resonance. The results show that the effect of parametric resonance is the largest 
at the half of the frequency of the oscillating boundary (u% = fi/2). This can be understood 
by considering the Mathieu equation 

x + /i 2 (l + ecosfit)x = 0, (5.1) 

where the parametric resonance takes place most strongly for Q = 2fi. In addition, in the 
case 7 > 2, we see the other resonance effects in addition to u>k = which is due to the 
effect of couplings with other mode frequencies in the cavity. 

We used the Floquet theory to examine the long-time behavior of the solutions and 
introduced the general scheme to find the solution that is available in the long time. For the 
simplest case, we found the characteristic exponents and the corresponding eigenvectors. In 
the parametric system, the stability-unstability structure is important because the resonance 
condition, uj n = Q — in ( |3.16| ), is hardly satisfied exactly in the experimental situation. 
In fact the condition of parametric resonance admits some discrepancy as seen from the 
solutions of the Mathieu equation. Therefore it is expected that the above simplest case is a 
good model to study the stability-unstability structure in the system of coupled parametric 
oscillators. We hope to report on this structure in a future paper. Finally we would like 
to mention that it is remained to solve the three term recurrence relation ( |4.13| ) in a future 
study. 
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